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Abstract 

The very high multiplicity (VHM) processes are considered to in- 
vestigate consequences of the dual resonance models reach (exponen- 
tial) mass spectrum. The virial decomposition is developed for de- 
scription of the produced particles short-rang (resonance) correlations. 
It is shown that the reach mass spectrum is able to do more flat the 
VHM distribution: a n = 0(e~ n ) up to multiplicities n ~ n 2 , n is the 
mean multiplicity. But for n >> n 2 one should expect a n < 0(e~ n ). 
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O ■ 1 Introduction 
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O . The last two decades of field theory development is marked by considerable 

efforts to avoid the problem of color charge confinement formulating a closed 
hadron field theory. The remarkable attempt based on the string model, in 
its various realizations Q. But, in spite of remarkable success (in formal- 
ism especially) there is not an experimentally measurable predictions of this 
^ . approach till now, e.g. [[|. 

The string model is a natural consequence of the old dual resonance model 
[0] and we hope that our toy approach includes main characteristic features 
of this model. We would like to describe in this paper production of 'stable' 
hadrons through decay of resonances. This channel was considered firstly in 
the papers ||. 

Our consideration will use following assumptions. 

A. The string interpretation of the dual-resonance model bring to the 
observation that the mass spectrum of resonances, i.e. the total number 
p(m) of mass m resonance excitations, grows exponentially: 

p(m) = (m/m ) 7 e /3 ° m , /3 = const, m > m . (1.1) 

Note also that the same hadron mass spectrum ( |1 . 1|) was predicted in the 
'bootstrap' approach ||. It predicts that 

7 = "5/2. (1.2) 
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Another assumptions are based on the ordinary (resonance «-> Regge pole) 
duality. 

B. The 'probability' of mass m resonance creation a R (m) has the Regge 
pole asymptotics: 

a R (m) = g R — , m>m ^ 0.2Gev, g R = const. (1.3) 

It was assumed here that the intercept of Regge pole trajectory a R = 1/2. 
So, the meson resonances only would be taken into account. 

C. If a R (m) describes decay of mass m resonance on the n hadrons, then 
the mean multiplicity of hadrons 

n R (m) = (L4) 



Following to the Regge model, 



m 
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n R (m) = n^ln^. (1.5) 

D. We will assume that there is a definite vicinity of n R {m) where <r R (m) 
is defined by n R {m) only. I.e. in this vicinity 

a R (m) = a R {m)e- nR{m \n R {m)) n /n\. (1.6) 

This is the direct consequence of the Regge pole model, if m/m is high 
enough. 

The connection between S-matrix approach and the real-time statistics 
(finite temperature field theories) f| will be used to formulate our model 
quantitatively. This interpretation will be useful since it allows to formu- 
late the description in terms of a few parameters only. All this statistical 
parameters are expressed through created particles energies and momenta. 

We will use also the virial decomposition technique. It was extremely 
effective for description of the phase transitions critical region, where the 
correlation radii tends to infinity. The Mayer' decomposition over 'connected 
groups' is well known in this connection |5]]. 

Following to our idea, we will distinguish the short-range correlations 
among hadrons and the long-range correlations among resonances. The 'con- 
nected groups' would be described by resonances (strings) and the interac- 
tions among them should be described introducing for this purpose the cor- 
relation functions among strings. So, we will consider the 'two-level' model 
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of hadrons creation: first level describes the short-range correlation among 
hadrons and the second level is connected to the correlations among strings. 

2 Asymptotic estimations 

Our purpose is to investigate the role of exponential spectrum (|1.1|) in the 
asymptotics over multiplicity n. In this case one can valid heavy resonances 
creation and such formulation of the problem have definite advantage. 

(i) If creation of heavy resonances at n — > oo is expected, then one can ne- 
glect the dependence on resonances momentum qi . So, the 'low-temperature' 
expansion is valid in the VHM region. 

(ii) Having the big parameter n one can construct the perturbations ex- 
panding over 1/n. We will see that there is a wide domain for n, where one 
can neglect resonance correlations. 

(iii) We will be able to show at the end the range of applicability of this 
assumptions. For this purpose following formal phenomena will be used. Let 
us introduce the 'grand partition function' 

E(z, s) = ^2z n a n (s), E(l,s) = cr tot {s), n < x/i/m = n max (s), (2.7) 

n 

and let us assume that just H is known. Then, using the inverse Mellin 
transformation, 




This integral will be computed expanding it in vicinity of solution z c > of 
equation: 

d 

n = z— lnS(z,s). (2.9) 

It is assumed, and this should be confirmed at the end, that the fluctuations 
in vicinity of z c are Gaussian. 

It is natural at first glance to consider z c = z c (n, s) as the increasing 
function of n. Indeed, this immediately follows from positivity of a n (s) and 
finiteness of n max (s) at finite s. But one can consider the limit mo — > 0. 
Theoretically this limit is rightful because of PCAC hypotheses and nothing 
should be happen if the pion mass mo — > 0. In this sense H(z, s) may be 
considered as the whole function of z. Then, z c = z c (n, s) would be increasing 
function of n if and only if S(z, s) is regular function at z = 1. 
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The prove of this statement is as follows. We should conclude, as follows 
from eq. (|2.9|) , that 

z c (n, s) — ► z s at n — ► oo, and at s — const, (2-10) 

i.e. the singularity points z s attracts z c in asymptotics over n. 

If z s = 1, then (z c — z s ) — > +0, when n tends to infinity [|7|]. But if z s > 1, 
then (z c — z s ) — > — in VHM region. 

On may find the estimation: 

_ 1 ^ ^(s) = ^ + 
n atot(s) 



where z c is the smallest solution of ( |2.9|) . It should be underlined that this 
estimation is independent on the character of singularity, i.e. the position z s 
only is important. 



3 Partition function 



Introducing the 'grand partition function' 
means that 

E(z,/3) 



1 



?s?/n 



In 

&tot(S 

'•' ( d 3 qidm^{qi, z)e~ f3e ^ 



the 'two-level' description 

-.-/3F{z,s) = 
N k { qi ,q 2 ,...,q k ;P), (3.12) 



i=l 



(2tt) 3 2^ 

where = yjqf + mf. This is our virial decomposition. Indeed, by definition 

^(z, s ) l^=i = vtot(s). (3.13) 
The quantity £(q, z) may be considered as the local activity. So, 

5E 



cr tot Ni(q) 



(3.14) 



5=1 



So, if decay of resonances form a group with 4- momentum q, then Ni(q) is 
the mean number of such groups. The second derivative gives: 

5 2 E 

~ Vtat{N 2 (qi, q2) - Ni(qi)Ni(q 2 )} = o tot K 2 {q x ,q 2 ) 

(3.15) 



^(qi, z)Si(q 2 , z) 
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where K 2 (qi,q 2 ) is two groups correlation function, and so on. 

The Lagrange multiplier (3 was introduced in ( 3.1 2| ) to each resonance: 
the Bolzmann exponent exp{— f3e} takes into account the energy conservation 
law J2i e i — E, where E is the total energy of colliding particles, 2E = y 7 ? in 
the CM frame. This conservation law means that (3 is defined by equation: 

yi=^lnH(z,/3). (3.16) 



So, to define the state one should solve two equations of state (|2.9|) and 

The solution f3 c of the eq.( |3.16| ) have meaning of inverse temperature of 
gas of strings if and only if the fluctuations in vicinity of (3 C are Gaussian. 

On the second level we should describe the resonances decay onto hadrons. 
Using ( |1.6|) we can write in some vicinity of z = 1: 

Z( q ,z) = Y,z n *N(q)=g R (—)e {z - lMm \ m=\q\. (3.17) 

The assumptions B and D was used here. 
So, 

- $T(z, s) = ^2 Ylidm^lrrii, z)}N k (m 1 , m 2 , m k ; (3), (3.18) 

k t=l 

where £ was defined in (|3.17|) and 

N k (m 1 ,m 2 , ...,m k ;P) = J J[l — 2e (q ) — J Nk ^ qii ?2 ' *"' qk] mi ' m2 ' •••> m ^)- 

(3.19) 

Assuming now that \qi\ « m are essential, 
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N k { mi ,m 2 , m k - (3) ~ N' k ( mi ,m 2 , ...,m k ) J[ U^e'^ J (3.20) 
Following to the duality assumption one may assume that 

k 

N' k ( mi ,m 2 ,...,m k ) = N k (m u m 2 ,...,m k )Y[{m]e^} (3.21) 

i=l 
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and iVfc(mi, 777.2, ■ ■■,' m k) is slowly varying function: 

N k (m 1 ,m 2 , ...,m k ) ~ C fc 
In result the low-temperature expansion looks as follows: 

- (3T{z, s) = £ 2k/2m ^ k 9 / R / Ck { j°° dmm^e^W-^J* . 

(3.22) 

We should assume that (/3 — /3 ) > 0. In this sense one may consider l//3 
as the limiting temperature and above mentioned constraint means that the 
string energies should be high enough. 



4 Thermo dynamical parameters 

Remembering that the position of singularity over z is essential only, let us 
assume that the resonance interactions can not renormalize it. Then, living 
first term in the sum ( |3.22j ), 

_ pjr^ s) = m °£^ f°° dm{m/m Q y^' 2 e^ nR ^-^-^ m . (4.23) 



^3/2 
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We expect that this assumption is hold if 

nmn n 

n -» oo, s oo, —=- = << 1. (4.24) 

So, we would solve our equations of state with following 'free energy': 

-f&M = ^ f^)(™)^e-^), (4.25) 
where, using (|1.2|), 

7' = 7 + 2(z - 1)%q + 5/2 = 2{z - l)fi* A = m (/3 - /3 ) > 0, a = const. 

(4.26) 

We have in terms of this new variables following equation for z, 

_ 2^^1(7^) 
n - z pz/2 dl , A y • ^0 
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The equation for (3 takes the form: 



_ amor(V + l,A) 

"-max - ^ 3/ / 2 ^ 7 '+l ' 

where n max = {y/s/mo) and r(A, 7') is the incomplete T-function: 

r°° 

r( 7 ',A) = / dxx"i- l e- x . 

J A 

5 Asymptotic solutions 

Following to physical intuition one should expect the cooling of the system 
when n — > 00 (at fixed y/s) and heating when n max — > 00 (at fixed n). But, 
as was mentioned above, since the solution of eq. (|4.28| ) (3 C is defined the value 
of total energy, one should expect that (3 C decrease in both cases. So, the 
solution 

dA dA 
A c > 0, — - < at n -> 00, — - < at s -> 00 (5.29) 
an as 

is natural for our consideration. 

The physical meaning of z is activity. It defines at (5 — const the work 

needed for one particles creation. Then, if the system is stable and H(z, s) 

may be singular at z > 1 only, 

dz c dz c 

^— > at n — > 00, ^— < at s — > 00. (5.30) 
an as 



One should assume solving equations Q4.271) and (|4.28|) they 



^ +1 |,^^<<r( 7 : + l,A c ). (5.31) 

This condition stands the physical requirement that n « n max . In opposite 
case the finiteness of the phase space for m 7^ should be taken into account. 

As was mentioned above the singularity z s attracts z c at n — > 00. By this 
reason one may consider following solutions. 

A. z s = 00: z c » A, A << 1. 

In this case 

A- 7 'r( 7 ',A) ~e 7 ' ln(7 ' /A) . (5.32) 
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This estimation gives following equations: 

n = Ci7'ln( 7 7A)e 7 ' ln(7 ' /A) , 

= 2 A 7 'ln(^)«l, (5.33) 

where Cj = 0(1) are the unimportant constants. The inequality is conse- 
quence of ( [5.311 ). 

This equations have following solutions: 

A c ~ << 1, i c ~ Inn >> 1. (5.34) 

Using this solution one can see from ( |2.11| ) that it gives 

a n < 0(e" n ). (5.35) 
B. z s = +1: z c — > 1, A c << I. 

One should estimate 1^(7', A) near the singularity at z = 1 and in vicinity 
of A = to consider the consequence of this solution. Expanding 1X7', A) 
over A at 7' — > 0, 

IYY, A) = r(V) - A 7 V A + 0(A y+1 ) ~ - + 0(1). (5.36) 

7 

This gives following equations for 7': 

w= yin(l/A)-l eVln(1/A) _ {537) 

i 

The equation for A has the form: 

n max = C^'+W/A). (5.38) 

Where C[ = 0(1) are unimportant constants. 
At 

< 7'ln(l/A) - 1 « 1, i.e. at ln(l/A) << n « ln 2 (l/A), (5.39) 
we find: 

^Kikr < 5 - 40 > 
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Inserting this solution into ( [5.38 



A c ~ — . (5.41) 

It is remarkable that A c in the leading approximation is n independent. By 
this reason 7^ becomes n independent also: 

^ 5*b ■ z -= 1 + nd^y (5 ' 42) 

This means that 

a n = 0(e- n ) (5.43) 

and obey the KNO scaling with mean multiplicity n = \n(n max ) . 



6 Conclusion 

Comparing A and B solutions we can see the change of attraction points with 
rising n: at n ~ n 2 (s) = n^ln(\/s/mo) the transition from (|5.43 ) asymptotics 



to ( |5.35| ) should be seen. At the same time one should see the strong KNO 
scalings violation at the tail of multiplicity distribution. 

The comparison of considered above model with experimental data at 
moderate and high energies will be given in subsequent papers. 

We have neglect the strings interactions and the final state particles inter- 
actions deriving this results. This assumption seems natural since z c — 1 << 1 
is essential at n << n « n 2 . By this reason one can neglect higher powers 
of (z c — 1) in expansion of ln£(z, to) over (z c — 1). Therefore, describing 
£(x,m) we may restrict ourselves by Poisson distribution (|I~6|). 

At the same time, at first glance, we can not neglect in ( p.22| ) the con- 



tributions with k > 1 in the moderate region n « n « n 2 . Indeed, fc-th 
order in (|3.22|) is 



~ Y\i c , A c ) ~ (l) fc ~ (lnA c ) fc ~ (\nn max ) k ~ (Hs/m 2 )) k » 1 

Nevertheless it can be shown that the higher terms with k > 1 can not change 
our semiqualitative conclusion. This question will be considered later. 
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